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Abstract. In this paper we consider the dynamics of dislocations with the same Burgers 
vector, contained in the same glide plane, and moving in a material with periodic obstacles. We 
study two cases: i) the particular case of parallel straight dislocations and ii) the general case 
of curved dislocations. In each case, we perform rigorously the homogenization of the dynamics 
and predict the corresponding effective macroscopic elasto-visco-plastic flow rule. 



1. Introduction 

In the recent years, an important effort has been done, both to improve the methods to compute 
discrete dislocation dynamics (see for instance the book of Bulatov and Cai [1] and the references 
therein) and also to connect them to continuum models of plasticity in crystalline solids (see 
for instance Fivel et al. [2] and more recently Hoc et al. [3]). Although continuum models 
of dislocations are known since the 50's (see Kroner [4, 5]), the dynamics has been taken into 
account only recently : see Groma et al. [6, 7] in 2D (and their mathematical studies in [8, 9]), 
Hochrainer et al. [10], and Monneau [11] in 3D. The goal of our work is to present, on a particular 
example, a rigorous justification of a continuum model with densities of dislocations bridging 
the gap with dislocation dynamics at the microscale. Indeed for a very special geometry, we 
are able to deduce by homogenization, the macroscopic elasto-visco-plastic flow rule relating the 
plastic strain velocity to the shear stress. The full technical details are presented in [12]. 



2. Homogenization of straight dislocations 

In this section, we consider the case of parallel straight edge dislocations with the same Burgers 
vector b = be x with b > 0, where (e x ,e y ,e z ) is an orthonormal basis with corresponding 
coordinates (x,y,z). All these dislocation lines are assumed to be contained in the same glide 
plane (x, y) and to move in this plane. 



2. 1 . The microscopic model for straight dislocations 

Because of our assumptions, for every integer i G Z, we can simply describe the position of the 
i-th dislocation by its real abscissa that we call Xi{t) where t is the time. We want to take into 
account the interactions of each dislocation with other defects in the crystal, that constitute 
obstacles to their motion. Those obstacles can be for instance other pinned dislocations or 
precipitates. In order to simplify the analysis, we will assume that these obstacles are periodically 
distributed, of spatial period A. In our model, those obstacles will be simply modeled by a smooth 
periodic potential V per satisfying V per (x + A) = V per (x). Then the energy of the system is the 



sum of two contributions: the interactions of each dislocations with the periodic potential and 
the sum of the two-body interactions between dislocations associated to a pair potential V. The 
energy of a set of dislocations is then given by 



E = ^ V per \ Xi ) + ^ V ( x i ~ x j) with V(x) = -j2bln\x\ and p, 



2vr(l - v) 

where the constants /x and v are respectively the shear modulus and the Poisson ratio. Remark 
that the force —V'(x) is then the usual Peach-Koehler force created at the point x by an edge 
dislocation positioned at the origin. 

We then consider the fully overdamped dynamics, where the velocity is proportional to the 
force, i.e. 

B-^ = -V Xi E + r ext (2.1) 

where B is the viscous drag coefficient and the force is on the right hand side. The first 
contribution to the force is a term deriving from the energy and r ex t is a real exterior applied 
shear stress, that can be seen as a driving force of the system. A natural question is then: what 
is the macroscopic behavior of this system ? 

In order to answer this question (which is done in Theorem 2.1), we have to introduce the 
plastic strain. To each dislocation is associated a three-dimensional displacement in the crystal, 
whose plastic strain is localized in the glide plane z = and is equal to 7<5o(^) where 5o is the 
Dirac mass. For instance, for a dislocation Xj, the intensity 7 (that we continue to call plastic 
strain) is equal to —bH(x — xi) where the Heaviside function H(x) is equal to 1 for positive x 
and zero otherwise. Here the sign defining the plastic strain is such that the quantity 7 increases 
when xi increases. Then the total plastic strain can be written as 



7 (M) = -b^H{x - Xi{t)). 



2.2. The normalization procedure 

We are now interested in the behavior of the system at a macroscopic scale A such that 
A >> A = Xb where A > 1 is a fixed ratio. Then we introduce several dimensionless quantities. 
We call x and t the normalized spatial and time coordinates at the macroscopic level, and 
introduce a parameter e and the associated normalized macroscopic plastic strain 7 e such that 
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with 7 £ (x, 0) = e 
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(2.2) 



where [•] is the floor function, 70 is a given function and BA/fi is a typical macroscopic time 
deduced from equation (2.1). Remark that e can be very small in our application (for instance 
£ ~ 1(T 6 if b ~ lCT 9 m and A ~ l(T 3 m). 

We expect that the macroscopic behavior of the model is well described by the limit 
macroscopic plastic strain 7°(x,t) o/7 e (x,t) as e goes to zero. 



2.3. Heuristics for the macroscopic stress field 

In this subsection, we want to give heuristic expressions of the normalized dislocation density 
and the macroscopic stress field, in terms of the limit macroscopic plastic strain. 

We remark that the gradient of the map x 1— > — 7 £ (x/A, t)/e is a sum of Dirac masses, 
and then the number of dislocations in a large segment of length Ax is formally given by 
— J Q Ax jK7$f( x /A-, t) dx. This shows at least formally that the dislocation density can be 



estimated as p(x,t) = —j\^^{ x ^t)- Then the total stress on the right hand side of (2.1) 
can be formally described at the macroscopic scale by 



r+°° dx' ^7° 

T = T ext + T sc with T sc (x,t) = -fi _ _ _ / -^—(x',t) (2.3) 

J X X (jX 

where we take the principal value in the integral defining the self-consistent field t sc . This 
expression can be deduced from the equation T sc (x,t) = —(V * x p)(x,t), where * x denotes the 
convolution with respect to the variable x. Remark also that the expression (2.3) of r sc is known 
to be the resolved shear stress created by the normalized dislocation density 

where for instance p° = 1/A when there is one dislocation by spatial period A. In particular, we 
see that r sc keeps the memory of the long range interactions between dislocations. 



2.4- The homogenization result 

We expect that the effective equation satisfied by the limit 7 can be written 
( <9t° 

— U = /(p°,r), for all xeK, te(0,+oo), 
< dt (2.5) 

k 7°(x,0) = 7o(x) for all x£R 

where p° is given in (2.4) and r in (2.3). Then our main result is: 
Theorem 2.1 (Homogenization of straight dislocations) 

Assume that the initial data 70 is non- decreasing and satisfies I70I + |7ol + l7ol — C for 
some constant C . Then for any C 2 periodic potential V per , there exists a continuous function 
f : M? — > M. such that r 1— > f(p ,r) is nondecreasing. And there exists a unique viscosity solution 
7 of the equation (2.5). 

Moreover, under the assumptions and notation of this section, there exists a unique solution 7 s 
associated to the dynamics (2.1) with initial data given in (2.2), and j £ converges to 7 locally 
uniformly on R x [0, +00). 

This result is proven rigorously in [12] in the mathematical framework of viscosity solutions 
(see for instance Crandall, Ishii, Lions [13] for an introduction to this theory). We explain in 
the next section how we compute the function /, which keeps the memory of the short range 
interactions between the dislocations and the periodic potential V per . 



2.5. Computation of f using Orowan's law 

In this subsection, we briefly explain (without any justifications) how to compute the function 
/. We refer the reader to [12] for the proofs of those results. 
Case A: V per = 0. 

In this special case, we can show that 

f(p jT )=p° v with v= T - (2.6) 

A* 

which is nothing else than the normalized Orowan's law giving, in a dimensionless form, the 
plastic strain velocity as the product of the normalized dislocation density p° and the normal- 
ized mean velocity v of the dislocations. 



Case B: General periodic potential V per . 

In that case, the function / can be computed using the following two steps. 
Step 1. 

For i £ Z, we look for solutions to (2.1) of the following special form 

( vt i \ - - 

x,i{t) = b ■ h (— + — J , with h(a + A) = A + h(a) for all a€R 

for some constant v and for a function h which is called a hull function. Both v and h have 
to be determined. Because of the convexity of the two-body potential V outside the origin, it 
is possible to show that the constant v exists and is unique. Moreover this constant v can be 
interpreted as the mean velocity of each dislocation. 
Step 2 . 

We simply define f(p°,T ex t) using the normalized Orowan's law as in (2.6), but with the 
normalized velocity v replaced by the constant v = M-v. 



2. 6. Numerical computation of f 

We present numerical simulations for the computation of the function /. We work with 
dimensionless quantities: X = l = X = b = B = p. We put initially N dislocations in an 
interval of length I = 10 which is repeated periodically. Therefore this interval contains I times 
the period of the periodic potential that we choose equal to V per (x) = sin(27rx) with A = 3. 
We discretize the ODE system (2.1), using an explicit Euler scheme with a time step At = 0.01. 
We compute numerically the mean velocity v of the dislocations after a final time T = 1000. We 
then set / = p°v with p° = N/l. We do the computation with N = 1, 200 and < r ex t < 9 
with Ar ex t = 2§g. Remark that we can restrict our computation for positive T ex t, because we 
have f(p°, —T ex t) = —f(p°,T ex t), from the symmetry of the potential V per in our problem. The 
level sets of the function / are represented on Figure 1. In order to have a better view of the 
set where / = 0, this set is conventionally represented in Figure 1 with artificial negative values 
of /. We remark that this figure shows in particular a collective behavior of the dislocations: 
higher is the density of dislocations, then easier the dislocations move above the obstacles. 

Figure 2 shows the map r ex t i— » f(p°,T ex t) for p° = N/l with N = 1,10,20. We see in 
particular that for r ext under a threshold (that depends on the dislocation density p°) the 
function / vanishes. 
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Figure 1. Level sets of the effective Figure 2. For N = 1,10,20, graph of the 

f {N/l, T ex t) with iV on abscissas and r ex t on map r ext h- > f(N/l,T ex t) 

ordinates 



3. Homogenization of curved dislocations 

In this section, we very briefly generalize the previous analysis to the case of curved dislocations 
all contained in the same plane (x, y) with the same Burgers vector b = be x with b > 0. 



3. 1 . The microscopic model for curved dislocations 

For ieZ, the motion of the i-th dislocation curve Ti(t) at the point X E M 2 is given by its 
normal velocity V defined by 

B-V(X,t)=r^(X) + ^F 3 (X,t) (3.1) 

3 

where Fj(X,t) is the resolved Peach-Koehler force created by the dislocation Tj(t) at the point 
X. Here r per is a smooth periodic function satisfying r per (X + Xk) = T per (X) for all k G 1?, 
which represents the periodic obstacles to the motion of the dislocations and can also include 
the exterior applied stress. To give the expression of this force, it is convenient to introduce a 
continuous function j(X, t) such that each dislocation curve Tj(t) can be seen as the level set 
l(X, t) = jb (when this level set is non-degenerated). Then a good approximation is given by 

Fj(X, t) = \f dZ J(X - Z) sign( 7 (Z, t) - jb) 

where, in the integral, the sign function takes values —1,0, 1. Here the kernel J is smooth and 
satisfies for a cut-off radius R = Rb with R > 1 fixed: 



J(-X)=J(X)>0, and J(X) = J^X) := -^—g ( ^- ) for \X\ > R > 



where for isotropic elasticity with X = (x,y), we have g = ^ { ~ ^ W+y 1 ^ ^ } w ^ n 

(3 = Remark that this formula allows to describe with the same formalism edge, screw and 
mixed dislocations (see for instance [14]). We also define the plastic strain 7 as 



7 = 6 



where we recall that [•] is the floor function. Then we proceed as in the previous section and 
define 



(3.2) 



3.2. The homogenization result 

We expect that the effective equation satisfied by the limit 7 of j £ can be written 
-^ = /(-V 7 °,r sc ), for all XeR 2 , te(0,+oo), 
, 7 °(X, 0) = 7 (X) for all iGl 2 

with 

t sc (X, t)= [ dZ J^{X - Z) 7 °(Z,t) 

where we take the principal value of the integral. Remark that this expression of t sc is consistent 
with the one given in (2.3) in the special case where j°(x,y,t) is independent of y. Then we 
have 



(3.3) 



Theorem 3.1 (Homogenization of curved dislocations) 

Assume that the initial data satisfies 70 + IV70I + \D 2 jq\ < C for some constant C. Then 
for any C 2 periodic function r per , there exists a continuous function f : M 2 x R — > R suc/i £/iai 
T l— ► /(') r ) nondecreasing. And there exists a unique viscosity solution 7 0/ i/ie equation 
(3.3). 

Moreover, under the assumptions and notation of this section, there exists a unique solution 7 s 
associated to the dynamics (3.1) with initial data given in (3.2), and 7 s converges to 7 locally 
uniformly on M 2 x [0, +00). 

4. Conclusion 

The main result of our work is the justification of the elasto-visco-plastic flow rule by the 
homogenization of the dynamics of dislocations with the same Burgers vector, moving in the 
same glide plane with periodic obstacles. Even if this geometry is very particular, this is, up 
to our knowledge, the first rigorous result in this direction. We also explained how to compute 
the flow rule, and presented numerical results. The proof of the homogenization for straight 
dislocations uses strongly the local convexity of the two-body potential V (which is equivalent 
to the non- negativity of the kernel J in the case of curved dislocations). 

Remark that for the same dynamics, it is possible to find non-convex potentials V, for which 
there is no homogenization. For a general geometry, there is in general no hope to find any 
convexity argument to justify homogenization. On the contrary, it seems reasonable to think 
that homogenization could arise in general, if we assume moreover that the dynamics is modified 
by the addition of a small random noise. But this is still an open problem to investigate. 
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